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MODEL TASK FOR THE DYNAMICS OP AN UNDERWATER TWO-LEGGED WALKER 
V. V. Beletskly, V. V. Golubkov, Ye. A. Stepanpya 

Academy of Sciences USSR Institute of Applied 
Mathematics imenl M. V. Keldysh 

Kecently Investigations of the dynamics of walking of anthro- /h* 
pomorphlc apparatuses and man have been strongly developing. This 
direction Is connected with the maturing possibilities of creating 
Integral robots and anthropomorphic autonomic systems of a type of 
exoskeletons and pressure suits for operation In experimental conditions. 

The proposed work borders on the Indicated Investigations. A 
model task of two-legged underwater walking Is examined. This task 
allows for the establishment of several characteristics of the walking 
of two-legged underwater apparatuses or pressure suits. The under- 
water walking device Is represented by a substantial sphere, which 
moves on dual-member rigid legs under the action of momentums which 
are located In the Joints of the legs. The legs of the apparatus are 
round cylinders, the joints are hinged between themselves and with 
the sphere. 

The dynamics of this system are Investigated with the calculation 
of the buoyancy of Archimedes, which acts on the sphere, and the force 
of hydrodynamic resistance, which acts on the sphere and legs. For 
the assigned walking velocity of the apparatus, the compensating 
vibration of the housing, momentums of force In the Joints and the 
key reactions are determined. Several questions of stability are 
investigated. A comparison of underwater and terrestrial walking Is 
given. 

1 . Definition of the Task. The Dynamic Model of a Walking Apparatus 
and Equalization of its Motion. 

Plane-paraiiel motion of a walking apparatus In a liquid is 
•Numbers in the margin indicate pagination in the original text. 
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examined. The apparatus consists of a rigid housing In the form of 
a sphere and a pair of legs which are mounted to the housing at one 
point of support II (Pig. 1). 

For the description of this movement, let us introduce the 
stationary, right quadrate Descarte system of coordinates Nxyz, the 
plane Nyz which coincides with the plane of motion of the apparatus. 

Let us make the following assumptions: 

1. The legs of the apparatus are Identical and each of them 
consists of two joints — a thigh and shank. The thigh and shank are 
round cylinders, which have corresponding lengths a and b and 
diameters d^ and dg . 

The legs are weightless and sufficiently thin, l.e. 

where m^ and m, are the masses of the thigh and shank. 

3 . The attachment of the shank to the thigh, and the thigh to 
the body (housing) Is jointed with one degree of latitude. 

. The plane motion of the apparatus takes place In such a manner 
that the center of the sphere 0, the point of support n and the legs 
are found In the plane Nyz. 

5 . The distribution of mass on the sphere and inside It is such 
that the plane Nyz Is the plane of dynamic symmetry of the sphere 
(apparatus), and the center of mass C of the apparatus lies on a line 
which passes through the center of the sphere and the point of 
support . 

6. The force of gravity Is directed along z. 

7 . The forces of friction In the hinges (joints), the lifting 
forces of the leg joints and forces of viscid friction, which 


act on the apparatus from the side of the liquid, are negligible. 

8. The surface with which the legs of the apparatus contact is 
absolutely uneven. 

9. The control by the action of the apparatus is realized with 
the aid of guiding momentums affixed in its Joints. 

10. Acceleration and the speed of any point of the apparatus 
are finite. 

With the assumptuions made, let us determine the task of the 
motion of the walking apparatus in a liquid in the following 
manner. Let 

1) the motion of the point of support be known and take place 
uniformly and rectlllnearly along axis y (comfortable motion [1]): 

2) the legs of the apparatus move along a given regular marked 
path, placed in the plane Nxy on the line Ny (In other words, the 
apparatus walks along plane Nxy along the line Ny with a constant 
pace): 

3) the step of the apparatus is single-supported, periodic and 
repetitious [1,2]: at each point of time the apparatus is supporteJ only 
on one leg; the movement of the legs of the apparatus take place 
periodically with a period of 2T (T is the length of one step); and one 
leg exactly repeats the motion of the other with a delay of T; 

i<) the motion of the end of the moving leg is specified in the 
form of a clear function of time. 

The position of the body of the apparatus in space will be 
characterized by two Descartes points of support yj^, Zj^ and angle y» 




3 



Figure 1 


and the position of the legs 
by the angles a^, o^, 

(Pig. 1). 


Let us consider that the 
contact of each foot with the 
surface is a point and all the 
action of the surface on the 
foot is reduced to a single ^ 
force (the force of reaction) 

(i ■ y, ff), applied to 

the foot at the point of 
support . 

From here on, through 
u/"), ana u/">. 

we will designate the momentums 
of the force of control, which 
acts from the side of the body 
on the thigh and from the side 


of the thigh on the knee respectively In the support and moving legs. 

By the establishment of the task, the coordinates of the point of support 
yjj, Zjj and the angles a^, B^, a^, B^, which define the position of the 
leg in space, it is known to us as a function of time. 

The forces of reaction and also of control u^ appear as unknowns. 


Applying to the walking apparatus the principle of the elimination 
of connections and using the theories concerning the change of the 
quality of the motion of the apparatus and its kinetic momentum relative 
to the point of support, we receive the equation for the definition of 
unknown forces of reaction and the angle y. These equations are the 
actual development of the equation received in [1] and can be 






recorded In the form 


• • 

Ajt •{ wA* ^A A 


( 1 . 1 ) 


Here and henceforth the point denotes the dlf ferentlalizatlon 

(A) 

in time. In (1.1) is Indicated: I _ is the tesor of inertia of the 

-► 

apparatus at the point of support II: p , p are the radius-vectors 
leading from the point of support to the center o*' the sphere to the 
center of mass of the apparatus (they are, on the strength of assump- 
tion 5 on page 4, colinear); ?jj are the radius-veotors /9 

which lead from the beginning of the selected system of coordinates 
respectively to the point of support and to the supporting points of 
the feet; is the mass of the apparatus equal, on the strength 
of the weightlessness of the legs, to the mass of the body; Mj^p is 
the combined mass of the body-sphere; w is the angular velocity of 
the body; is the velocity of the point of support; is the radius 
of the sphere: P is the density of the liquid; and f'g are vectors 
of the buyolng force and force of the weight which acts on the body; 

^v, ^ the forces of hydrodynamic resistance of the liquid 

which acts respectively on the body (sphere) and the legs of the 
apparatus; i^., M are the momentums of force of the hydrodynamic 

resistance relative to the point of support, which act respectively on 
the body and the legs of the apparatus. 


In deriving equation (1.1) in a known approximation we 
considered the Instability of the streamlining of the body-sphere of 
the walking apparatus in its motion in a liquid by means of the intro 
ductlon of its combined masses [3] (the combined masses of the legs 
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were assumed to equal zero, which on the strength of supposition 
2 on page 4 is fully feasible). 

From here on only the single support motion Is examined. It Is 
convenient to introduce the Index "v" for all values which pertain 
to the support to the leg and the Index ”ir" for all values which 
relate to the moving leg (Pig. 2,3). Then on the strength of the 
establishment of the task and suppositions which have been made, we 
have (Fig. 1, 2, 3): 




, ( 1 . 2 ) 

it.o. 

Hmconst , • const , /f/ • const . 


In addition 

T-constJ>const,tj,-(f^y‘*tJ: 

if 
if 


A- 


\p^\*0. 


/II 


(1.3) 
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Pg are the absolute values of the vectors 






^Is the point of Inertia of the apparatus relative to the axis 
which passes through the point of support and parallel coordinates 
of the axis x; H Is the height of the point of support over the sur- 
face Nxy, hj^ Is the projection of the radius of the vector « j| on 


the axis 


y at the Initial moment of time t^: 

0 V 


of the point of support; « 


is the radius-vector 
Is the radius-vector of the end of the 


moving foot; h Is the value of the step of the apparatus; T Is the 
earlier Introduced time In which the apparatus makes one step; x 


±0 


•*0 


z are respectively the unit vectors of the coordinate axes x, y, 
z; [t-to/T] Is the entire part of the unit t-t_/T. 


With the written formula (1.2) it was assumed that at the initial 
moment of time one leg Is found at the beginning of the system of 
coordinates N and the apparatus has ceased to be supported on this 
leg (the leg has started moving), and the other leg is found on the /12 
axis at a distance h from point N and the apparatus has Just 
begun to be supported on this leg (the leg becomes supporting) (Fig. 3). 




Figure 2 Figure 3 
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Substituting (1.2)-(1.3) into (1.1) and writing out the obtained 
vector of the equation on coordinate axes, we obtain the equation 
which describes the change of the angle y with timet and the expression 
for the projection of the forces of reaction. 

* ^finp i (t)cos^) i‘ 

* •* ^jr] . 

‘ ■ 'Sir “ ^J^„p . 

(M % 

, h,.h.h, , 


The four last formulas of (1.^) are written with the calculation of /13 
connections which exist between Vjj, T and h 

h . V„T. 

This connection is the consequence of the periodicity of motion of 
the walking apparatus. 

With the absence of hydrodynamic forces and momentums the equations 
of (1.4) are converted to earlier knowns [1]. 
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The first of the equations (1.4) has the structure 


d.*') 


In the task which Is examined In the present work, the right part of 
this equation depends on the clearly forward moving time t with a 
periodic pattern so that f( Y» t 2T) ■ f(v» Y» t). Let us desig- 
nate the value of the phase changes at the Initial moment t^ through 
Yq» Yq* Then the general solution of the equation (1,4') may be 
written In the form 


' ^ ff* */• » ^)t 


cft(r0j0.uj) .4. 

W ' «-e. 


(1.5) 


According to the establishment of the task y(t) must be a 
periodic function with a period of 2T. This Indicates that from the 
series of solution of (1.5) we must choose only those solutions which 
meet the following condition of periodicity: 

W.. 

,/f*. (1.6) 

Let us assume that the conditions of periodicity of (1.6) may 
be satisfied by the choice of Initial givens y, y. Then the relation- 
ships of (1.6) must be examined as a system of two transcendental 
equations with two unknowns y^, y^. Finding the solution of the 
system of equations of (1.6), we find by that very solution the periodic 
solution of the differential equation In (1.4), and, as It has come A4 
to be, the periodic motion of the walking apparatus. Thus the task 
of finding the periodic solutions of the equation for y Is reduced 
to the solution of the angular task of (1.4), (1.6). Below will be 
shown that this angular task actually has the solution. 


2 . Finding the Control Momentums In the Joints 

For the realization of the Indicated above motion of the walking 
apparatus. It Is necessary to place corresponding momentums of control 
of the forces u^ , u^K) (1 ■ v , ir) In Its Joints. Let us find these 
momentums. According to the assumption, they are located only In the 
hip and knee Joints. Let us apply to the parts of the leg (from the 
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gears to the point of support) the principle of the freedom of 
movement of connections and the theory concerning the change of 
kinetic momentum (respectively relative to the knee and the point 
of support). 


Let 1 ■ V ,ir Is the Index respectively of the support and the 
moving legs. Let us take the 1-leg and having eliminated the body 
and the support surface, let us change the motion of the body on the 

thigh at the point n of support by the force of reaction and 

momentum u^ * relative to n and the motion of the surface on the 
foot of the force of reaction (Fig. ^). Applying the theory con- 
cerning the change of kinetic momentum of the leg relative to the 

point of support and using assumptions of the weightlessness of the 
leg and the flnlteness of velocity and acceleration of any of Its 
points, we obtain 

Here are the momentums of forces of hydrodynamic resistance, 

which act respectively on the support and moving legs. 


Considering that 


from (2.1) we find u^*‘^ and consecutively assuming 

1 ■ V and 1 ■ TT 


/lo 


( 2 . 2 ) 


Let us now apply the theory of kinetic momentum relative to the 

knee to the shank of the 1-leg. Substituting the action of the thigh 

4 (K^ (K) 

on the shank by the force of reaction R^^ and momentum uj , and 
the action of the support of the surface — of the reactions at the 
point of support and taking Into c 'nslderatlon the weightlessness 
of the shank and the flnlteness of velocity and acceleration of Its 
points, we obtain (Fig. 5) 




10 



Figure Figure 5 

This gives, with 1 ■ v and 1 ■ n respectively 

^jr ' 

Here and are the momentums of forces of hydro- 

dynamic resistance which act respectively on the shank of the 
supporting and moving legs and are taken respective to the thigh; 

6 v Is the angle B for the support leg. From the given formulas with 
the absence of hydrodynamic momentums the formulas of work are obtained 
[1]. The obtained formulas for the control In the Joints (2.2, 2.3) 
hold for the single-support motion In a smooth motion of the walking 
apparatus In a liquid. From these formulas It is seen that the con- 
trol In the Joints of the moving leg In an examined case, despite 
Its weightlessness, are not equal to zero. In contrast to the situation 
when the apparatus moves In a nonresisting environment [1]. In order 
to achieve movement of the leg In a resisting envlronmtnt In the thigh 
and shank, it is necessary to place compensating controlling points. 

They are, according to the value, equal to the momentum of the forces 
of resistance which act respectively on the leg and the shank. After 
the solution of the angular task for y and definition of reactions, /17 
according to formulas (2.2), (2.3) the controls in the Joints are 
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clearly computed (If trajectories are assigned of the motion of the 
points of support and foot of the moving leg, and also the path 
being followed). 

3. Computation of Forces and Momentums of Forces of Hydrodynamic 
Resistance which Act on the Walking Apparatus 

For the calculation of angle y It Is necessary to know how 
to compute the forces and momentums of forces of the hydrodynamic 
resistance which act on the walking apparatus. Let us assume that 
the fluid is incompressible; the forces of viscid friction are 
negligibly small In comparison with the forces of hydrodynamic pressure; 
the calculation of the hydrodynamic forces and momentums which act 
on the various parts of the apparatus may be arrived at independently, 
without calculation of the interference; the motion of the apparatus 
is quaslstationary in a sense that for the calculation of the forces 
of hydrodynamic resistance It is possible to use the formulas which 
are correct in the case of stationary streamlining of the bodies 
(the unstationarlness of the streamlining in a known approximation 
has already been calculated by us by means of Introducing the combined 
masses of the body of the apparatus). 

Using the assumptions made, let us calculate first the forces 
and momenuums of hydrodynamic resistance which act on the body of 
the apparatus. With calculation of these forces and momentums let us 
consider the movement of the body-sphere as being set, and the velocity 
of this motion equal to the velocity of the center of the sphere. The 
hydrodynamic effects which are connected with the revolution of the 
oncoming flow of fluid on the sphere reduces '■•o a single force which 
is applied to the center of the sphere. This force Is given by the 
formula [4] 

P'v j? "• (3*1) 

where C .(He) is the coefficient of resistance which is dependent /l8 

C (p ' 
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the Reynolds number Re; p and v are respectively the density and 
coefficient of the kinematic viscosity of the fluid. 

The dependence of the coefficient on the Reynolds number 
Is shown In the following table 


Table 4 



O.I 

I 

10 

10^ 

Iq3 

10^ 

10® 

I(^ 


245 

28 

4.4 

1. 1 

0,46 

0.42 

0,49 

0.14 


In the work the motion of the apparatus Is numerically 

Investigated with - 1 and V -Im/sec. 

ce o 


With these conditions, as the value shows 


^ 0-0t?3 



the Reynolds number has a great value. Therefore the value of the 
coefficient In the formula (3.1) in our case must be derived 
for great Reynolds numbers. 

As has been noted, (3.1) was received as a result of the solution 
of the task of stationary streamlining of the sphere by a viscid. 
Incompressible fluid; the hydrodynamic effects of the unstatlonary 
streamlining have been disregarded. However, these effects add to 
the simple value, for example, on the basis of the task concerning 
the periodic oscillatory movement of the sphere in a viscid fluid 
which was examined in the work [6]. In this work is noted the 
following expression for the momentum of forces which act from 
the side of the fluid on the sphere which achieves a periodic 
oscillatory motion. 




Mm 




(») 
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where n is the viscosity, fi is the angular velocity of the oscillations 

of the sphere, u is the frequency of the change of the angular 

velocity, R is the radius of the sphere. The formula (•) is valid /19 

2 

in the hypothesis pR u >> n, which in the given case is realized. 

The values show that the value according to (•) does not exceed 
3% from the base momentum of hydrodynamic resistance which is cal- 
culated on the basis of formula (3*1) • These formulas lead to the 
following expression for the momentum of the forces Mj^ and Mj^p, 
which are defined by the motion of the forces of hydrodynamic resistance 
on the body of the apparatus 

Ka CK > 


Thus , 

(ay^ sifty*//cosy* )J 

■ f • '3.2') 


Let us calculate now the forces and momentums of the hydrodynamic 


resistance which act on the legs 

I 

I 



Figure 6 


of the apparatus. Let us take any 
leg and examine an arbitrary section 
of this leg (Fig. 6). An infinitely 
small element of this section of the 
leg d^ presents Itself as a round 
thin cylinder. 

With motion of the legs of the 
apparatus the angles of incidence 
are close to 90°. Considering that 


! 


Il4 




it is possible to disregard the influence of viscosity on the stream- 
lining of uhe cylinder under the angle of incidence, which is close 
to a right angle, we determine that hydrodynamic force acts along the 
normal to the axis of the cylinder. The force which acts on the 
element of the length of the leg dC, is described by the formula 
" ^pP(Vjj/2)d*dC, where Cp is the coefficient of any resistance 
of the cylinder at the angle of incidence a * 90 , Vjj is the component 
of velocity in the direction perpendicular to the axis of the cylinder, 
d is the diameter of the cross section of the section of the leg. 

In a vector form this may be recorded in the following fashion 


Here C is the distance from the beginning of the section (point /20 

H) to the element dC ; are the velocities respectively of the 

^ ^ ^ ^ 

element ci^ and the point H; (V,n), (Vjj,n) are the scalar products 
respectively of the vectors ^,n, and n. The coefficient of 
resistance depends on the Reynolds number. 


where V is the characteristic velocity of the motion of the element 
of the leg. In a characteristic situation with motion of the apparatus 
in water V = V ~ Im/sec with d ~ 0.09 mm we will have 


— - 0.09 

^ " 0.0 17 3 • tO'*> 


a 0.5 'fO^ 


In conformance to the results of work [5] in a sufficiently large 

area of values of the Reynolds number being examined it is possible 

to accept C ~ O.G. 

P 
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Considering that the streamlining of each element of the section 
move independantly , we receive the following expression for the full 
force which acts on the section 


m 


(3.^0 


By an analogous method are computed the momentums of the forces 


of hydrodynamic resistance relative to the point of support (Pig. 6). 


- - { 4 dph fVJCf . . 


(3.5) 


Substituting n, p and Vj^ from (3-3) into (3*^) and (3.5) we find 

- j Cfpd[f(j/^ siiMf-i, cosv)tmJ. 

(3.6) 

ij . K,,/. * v„S' . 


From (3.6) we have 

*M^X* ■ - ^h) ^ ^ 

/ *3 \ 

= m - i Cp pdf (Ay ^ $in(p *• Fcos y>) . 

Let us now apply formulas (3*6) and (3*7) for the calculation of 
forces and momentums of hydrodynamic resistance which act on the thigh 
and shanks of the walking apparatus. 
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For the thighs of the apparatus, assuming In formulas (3*6) and 

(3.7) 

* /•«/ . t^a , 

%’%‘V.yo . 

Sf^*o , g^.o 


we receive the following formulas fo the forces and momentums 

( 6 ) ( 6 ) ^ 

^1 * ^IP ' hydrodynamic resistance, which act respectively on the 

thigh of the moving (1 ■ ir) and support (1 ■ v) legs. /22 






( 6 ) 


.W I . . jii 


t •, (Si 

- J * f^cosofi ) 

(i) a jtt* , 

ft “ S *^9i 

(Si a jS) a*rjSi,jSi. off* , (S>* rt>* 

-fi f9, (9,‘ * 

V*- y, eosof. , (ff'm s^*6f. a , i ‘if, X . 


(3.8) 


The formulas for the forces and momentums of hydrodynamic 
resistance which act on the shanks of the apparatus we receive from 
(3.6) and (3.7) by means of substitution 


act’ll**, fflm^'eosof‘*i*siffa; , 
• ost/for^ y* - aeoscf^ i * . 



As a result the following formulas are received for the forces 
and momentums M^* ^ , M^p* ^ which act on the shank of the moving 
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(1 ■ tt) and support legs (l"v) 


^ •'S/ . 

*’J sinfi> * f/cos/ii ) . 


/at a W (t) 



With the aid of formulas (3.8), (3-9) the forces of hydrodynamic /^ 
resistance and their momentums which act on the legs of the apparatus 
are calculated conclusively in the following manners 

Fi ? 'Stii jSut Tfi) ^ct 

aa aa 

/fj> * • Mf * Mf * Mff * Mf( ^ 

Thus, the formulas received (3. D» (3*2), (3.2'), (3.8)-(3*10) 
which fully allow calculation of the forces and momentums of 
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hydrodynamic resistance which act on the walking apparatus during /24 
Its movement In a liquid. These formulas together with (l.ll) allow 
the calculation of angle y which describes the motion of the body of 
the apparatus relative to the point of support. 

^ . Numerical Method of Solving a Nonlinear Task of Walking 

Solution of the task of walking Is not reduced to solution of 
the angular task of (l.A), (1.6). The right parts of the differential 
equation (1.4) are described by the formulas (3.2), (3*2'), (3.8)- 
(3*10). The differential equation for y Is significantly nonlinear 
and complex. Therefore for the solution of the angular task (1.4), 

(1.6) It Is necessary to use numerical methods. In the given work 
the numerical. Iterative, generalized method of cords Is used. Descrip- 
tion of this method Is given In the Appendix. 

The value of the functions y(Yq» Yq» ^T) and YqCYo* 

to* tQ + ?T)for the dlffer<“nt vslues of y© Yo» which must be known 
In the process of solving the angular task (1.4), (1.6) by the general- 
ized method of cords, were received by a numerical Integration of the 
differential equation (1.4) by the method of Runge-Kutt with a constant 
step. Here the necessary values of the angles av(t),6v (t), ait(t), 
Bn(t), which are formed by the thighs and the shanks of the legs with 
the axis z, were calculated along the assigned trajectories of the 
point of support and foot of the moving leg with the help of the 
algorithm used In work [2]. 

Like any local algorithm, the generalized method of cords re- 
quires knowledge of the zero approximation. As a zero approximation 
It is possible to use either the solution of the task of walking which 
Is received in a linear approximation In the supposition of an 
Insignificant amount y{or y-ii)* Y andVor the values Yq * o(Yo 
Yq “ 0. The Indicated zero approximations are sufficiently good, in so 
far as calculations have shown, the oscillation of the body of the 
apparatus along y and y which are obtained as a result of solving the /25 
nonlinear task of walking is expressed in ^iianalytlc form through 
quadratures. The solution of this task is given below. 
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5. Solution of the Task of Walking In a Linear Approximation 


Let us assume that the equations of (1.4) have such solutions 
which have y» Y and V sufficiently small. Then, applying the usual 
procedure of linearization to the formula for y I n (1.4), we 
receive the equation describing the small oscillations of the body 
of the apparatus: 

Here 

Pm , q . . 

f J: V J/ (5.1) 

^0 UU ft) * ^ 

In the equation (5.1) the member mj^(t) results from the action 
of the forces of hydrodynamic resistance on the legs of the apparatus, 
and the member mj^(t) — the action of the buoying force, the force 
of weight and the forces of hydrodynamic resistance on the body of 
the apparatus. 

The general solution of the equation (5.1) may be obtained In 
the form 


* ch(cj,T)*e''*^( sh(a)^ (r-r))~c/i(a)f (r- T)))] , 
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rJ ^ 

■ ^ / ^ *^(<i)f9(t$t ck(Vf 9(f)) f 

• AJ- • 

7* ^ 

/ e*""'‘'W6j, wrw-f ^/)7,rt,*f;o'f , 


0 

,w w 




• • tit} /«! 


* 'c^ro^ r; rr~^)) - 




«V« /ifV' ' 


•cju 


1 ~h 


l<^,i-(ip/><f,0*^, 

^‘>-^^o^e/,-c/,.Jjtet, 


< 0 


Off* 


et . ^(P,*mv^ 

^--v: • 


/'-/ . 




jfckft 




Here Yq and y^ the values of Y(t).and Y(t) at the initial 
point of time t^. 

If P* - Qy < 0, then as is evident from formulas (5.2), uy is the 
virtual value. In this case y (t) is calculated according to those 
formulas of (5.2), only in them it is necessary to substitute uy for 
|uy| , ch for cos and sh for sin. In addition, from the formulas 
of (5.2) it follows that in the linear approximation the solution of 
Y (t) is presented in the fox*m of the sums of two components, one 
of which Yjj(t) results from the action of exterior forces on the body 
of the apparatus, and the second Yjj(t) by the action of exterior 
forces on Its legs. 


In turn, each component is the sum of two parts — the periodic 
(Yjj^"^(t), Yu^”^(t)) with a period of T and a nonperiodic or aperiodic 

(Yx^®^t), y^^Ht), 


In order that the solution of v (t) be periodic, it is necessary 
and sufficient that Yj^^^^(t) + ■ 0» which may be achieved 

by a selection of corresponding initial data Yq and y^. The initial 
data Yq and y^ which provides the periodicity of Y(t) is defined 
clearly and equals 


^ W (trt 

(t,) , 

<») g Off 

•[f*e'^^(-ch(cOgr)i.^ sACcj^r))] 


" . J e ^ sh«i)0<r-^)) ^ 
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(5.3) 
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The periodic solution has the form 


^ «;-4 ^4 , 


(") 


(n) 


(5.M 


where Mt) “ y^ Mt) an given by the formulas of (5.2). 


Analysis of the formulas of (5.2) shows that If qY>0, then the 
solution of Y (t' ' i'.vtnptotlcally approaches the periodic solution of 
Y^"!t). Thlf hij’.calor that In the case of qy > 0 the periodic /29 

oscillatory r .i?lon of ti:e body of the apparatus Is stable. If 
QY < 0 then '. hia motion Is unstable. 


Let us ncue that by ”stablllty" and ’"Instability" Is understood 
here the usual stability and instability according to Lyapunov of the 
periodic so'jjtlon of Y(t) of equation (5.1) relative to the perturbation 
of AYq, Ay^ of the Initial data Yq» y^. In tenns of the Initial 
dynamic task this means that the question lies In the conditional 
Invariability of the assigned forward motion of the point of support 
of the lers and the legs them? -Ives . The result of this Invariability 
of the t ■‘slgni forward motion (and motion of the legs) Is the 
Immobilization of the dyr^ndc parameters which enter Into equation 
(5.1). So, for •a/.Jlu, the parameters considered Invariable are 
V — the velocity of the forward motion of the point of support of 
the legs of the apparatus, H ~ the height of the point of support 
of the legs above the surface and so forth. 
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EttabllfhsMnt of the task of motion atablllty on the whole 
(relative to the entire aerlea of perturbatlona) leada to the taak 
of control — the taak of full walking atablllaatlon of the apparatua. 
In the preaent work the taak of full atablllaatlon la not examined. 


The deacrlbed effect of atablllty doea not have a place In the 
analogor^ altuatlon with the motion of the apparatua In nonrealatlng 
aurroundinga [1]. The noted feature of apparatua motion In a liquid 
la explained by the almultaneoua action on the body of the buoying 
force and the forcea of hydrodynamic realatance (In the abaence of 
a liquid, theae forcea are not preaent). 


Becauae for the examined caae of body poaltlon of the apparatua, 

§ 

cloae to the upper vertical, we have > 0, then the condlton 


j: 


( 5 . 5 ) 


Indlcatee that the reducing momentum of the buoying force of the body /30 
relative to the point of aupport of the lega la greater than the 
momentum of Ita force of gravity. Thia condition provldea the atabl- 
llty (In the aenae Indicated above) of the periodic motion of the 
body with the upper — relative to the point of aupport of the lega 
— poaltlon of the center of maaeea of the apparatua. The action on 
the body of the hydrodynamic forcea of realatance, which are dlaalpatlve, 
exceed the atablllty at the aymptotic. 


If the momentuma of the forc' ' ' .ich act on the lege are 
conaiderably lesa In comparlroi: with the momentuma of the forcea which 
act on the body of the apparatua, then It la then poaelble to aet 
m^(t) * 0 and for v(t) the following approximate value la received 

which In an explicit manner la expreaaed through the elementcry 
functlona. In thla caae, aa analyaia of the formulae for y(t) ahow. 
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in order for the procedure of linearization to be valid » It la 
necessary to require that 


\(P*-Pg)h I 






Cg(fi. *p)v: 

->?|A 




(5.6) 


In an analogous manner, In a linear approximation It Is possible 
to examine the motion In a case when the body of the apparatus Is 
close to the lower vertical position, l.e. receive and examine the 
solution of the equation for ir(l.^) In the area y * ^ * Y " 0, 
considering Ay ■ y-it* Ay and A*y with small values. In this case, the 
equation and formulas describing the change of Ay with time, are 
received from the equations and formulas for y when y Z 0 (formulas 
(5.1)-(5.M, If set In them 

The conditions for the applicability of the linear theory In /31 

the examined case have that form (5.6), only y^ must bc> substituted 
for AYq ■ Yq- If. The qualitative derivation about the character of 
body motion relative to the point of support of the legs In the area 

Y ■ IT, Y ■ 0 remains the same as In the case of y “ 0, y ■ 0. Only 
the conditions of the acymptotlc stability of the oscillations of angle 

Y change. In this case. It turns out that the oscillations will be 
asymptotically stable if 


and unstable If 




Here it Is relevant to note that according to the base equation 
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of motion (1.1), the momentum of the forces of control act on the 
apparatus, which —through the support reaction — rtepends on the 
angular acceleration of the apparatus. Thus the angular acceleration 
enters into the equation of motion not only through the change of 
kinetic momentum of the apparatus, but also directly — through the 
controls. In the examined planar task this is equivalent to the 
introduction into the examination of any effective "introduced" 
momentum of inertia — of the coefficient entering into the equation 
of motion by a multiple factor with the ang'«lar acceleration. The 
effective momentum of inertia, generally speaking, depends on the 
angular position of the apparatus, that is, it is not constant. It 

can even change sign. In the examined case of small oscillations 

§ 

the effective momentum of inertia Jjj is a constant value. But the 
sign of this value — dependent on the value of the parameters — may 
be positive or negative, which in turn, acts on the stability of the 
body motion of the apparatus in the sense earlier indicated according 
to the condition introduced above. And namely, in the case of oscil- 
lations of the body near the upper vertical position, the value of 

is positive: Jjj’* = Jj^ + MpjjpH >0. In this case, as the condition /32 
of stability was indicated (5 .5) reduces to the condition of posltlve- 
ness of the numerator in (5*5). 

With oscillation of the body near the lower vertical position 
we have J* = Jj^ + so that Jjj > 0 with Jj^ > MpjjpH and j|| < 0 

with Jjj < M Pj^pH. The condition of stability for such a case has the 
form of (5.7) and is fulfulled if the numerator and denominator in 
(5*7) have different signs. 

It is necessary to underline that the conditional stability 
defined above of the periodic oscillations of the body of the apparatus, 
a'i also all motion of the body, Implies calculation of the controlling 
momentum in the hip Joint of the apparatus. The great number of 
possible controlling momentums answer the assigned motion, being 
calculated by the formulas of (2.2) and are distinguished one from 
the other only through the various motions y(t) of the body. 














According to (1.1) and (2.2) with the assigned motion of the legs of 

the apparatus, the controlling momentum u^^^(t, y»Y»V) may be varied 

only on account of the variation of Ay, Ay, A*y* of the oscillation of 

the body which Is also calculated with examination of the stability 

of the periodic oscillation Y(t) of the body of the apparatus. The 

corresponding variation of control Is defined by the formula 

» 

where the sign "plus” answers the variation near the value y = 0, 
the sign "minus" — the variation near y = tt. 

In concluding this paragraph let us examine one significantly 
nonlinear effect. Important for the understanding of the mechanism 
of underwater walking. The periodic solution of equation (1.4) could 
be sought in the form of a Poure series, preliminarily having expanded 
to series the periodic coefficients of this equation. Realization of 
the procedure is hardly efficient in view of the great amount of 
computation; but the main, constant member of the periodic solution 
which gives a median position of the body, with a good approximation 
received from equation (1.4) by the establishment of y = Yq= const /33 
and with the retention of only the main, constant, members in the 
expansion of the coefficients leads to the formula 

Here < M > is the median value of the summated momentum of the forces 
of hydrodynamic resistance. 

The condition |slnY^i 1 , which is 

_ (Pa~ P»)( hj- &h)*< ^ y 

^aPc-p*Po ( 5 . 10 ) 


presents a necessary condition of the existence of the sought 
periodic, compensating motion of the body, which is the condition of 
the occurrence of the walking being examined. In the absence of hydro- 
dynamic resistance <M_> =0 and “ 0, and the condition (5.10) gives 

L D 
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->j3f*‘ih th^ , 


(5.11) 

(here It is accepted The inequalities of (5.11) define 

the series of possible configurations of the walking of the apparatus; 
these configurations with the assigned height H of the point of support 
above the surface and the assigned length h of the step, are clearly 
defined by the support section h 2 « 


The condition (5.10) may be recorded so: 

( 5 . 12 ) 

As it Is possible to understand from the formulas written out earlier 
for the momentum of the forces of resistance, we have < Mj> >0; 
accepting for the values p ~ p , we receive from (5.12) 


(5.13) 

This means that the diapason of possible values of the support section 
h 2 in underwater walking displaces to the side lesser values in com- /3^j 
parison with the diaposon in terreatrlal walking, ^rom Pig. 3 it is 
possible to understand that this corresponas, in turn, to those walking 
conf i.r'urations in which the point of support of the legs is shifted 
farther forward in motion, so that the body would be "placing its chest 
the wateri* moving somewhat forward relative to the legs. This fact 
is accounted for in numerical calculations where value h 2 = 0 is 
accepted. With h 2 ® 0 the point of support of the legs at the moment 
of change of the support leg is projected precisely on the point of 
support. Pres»=‘ntatlon of the result of calculation is given in the 
following parc-graph. 


6. Results of Calculations 


Let us rtescribe several results of calculations. The motion of 
a walking apparatus In water was examined. For the base variation » 
the variant with the following values of pareuneters was chosen; 


Radius of the sphere of the apparatus 

Distance from the point of support to the 
center of the sphere 

Distance from the point of support to the 
center of masses of the sphere 

Weight of the apparatus 

Buoying force 

Effective weight of the apparatus 

Velocity of the point of support of the legs 

Height of the point of support above the 
supporting surface 


R . ■ 0.9 m 

C9 

■ 0.6 m 

p = 0.3 m 
= 3180 kg 
Pg = 3130 kg 

V ‘ 50 kg 

Vjj = 0.5 m/sec 
H » 1.5 m 


Length of step 
Length of time of step 
Density of water 


h = 1 m 

T = h/Vj^ = 2 sec 

2 4 

p = 104 kg sec /m 


In addition. In the base variant, the configuration of motion Is 
defined by the fact that the point of support of the legs at the /35 
moment of change projects to the point of support (h^ * h, h^ * 0). 

In Fig. 7 results are shown of calculation of the base variant. 

In this figure behavior with time of angle y Is shown ~ the digression 
of the axis of the apparatus from the vertical, the components Fly, 

R 2 of the reaction of support, the controls uj" In the thighs and 
shanks of the supporting and moving legs. Hereafter these functions 
as a set will be called the "dynamic characteristics" of the apparatus. 
Behavior of the apparatus which follows after the results shown In 
Fig. 7» seemed to little resemble the behavior of an analogous 


29 








apparatus In the absence of resisting surroundings. 


Comparison of the given dynamic characteristics with the dynamic 
characteristics of [1] shows a qualitative and quantitative difference — 
In the behavior of the horizontal component of the reaction of 
support. For a terrestrial apparatus the function R (t) Is almost 
linear and changes sign in the middle of the step. For an under- 
water device this function is significantly not monotonlc and In 
almost all parts of motion is positive. The maximum value of it 30 kg 
with an effective weight of 50 kg (by the effective weight of the 
apparatus is understood the difference between its weight and the 
buoying force of the body). This is a factor of 10 greater than for 
an equivalent terrestrial apparatus. Such behavior of the horizontal 


component of th<" support reaction is 



Figure 7* Dynamic characteristics 
of the apparatus in its motion on 
one step. = 3180 kg, Pg = 3130 kg 

V„ = 0.5 m/sec, R . = 0.9 m, 
n c 9 

p =0.6m, p =0.3m. 

o ^ 


explained by the fact that it 
compensates the resistance of 
the environment. The vertical 
component of reaction is close 
to a constant value equal to the 
effective weight (50 kg). The 
controlling momentums in the 
Joints of the support leg are 
significantly nonmonotonic, in 
contrast to the terrestrial case. 
Controls in the thigh and shank 
of the support leg attain a max- 
imum (by the module) ~ kg m. 
Controls in the moving leg are 
much less (of an order of 5 kg m) ; 
they are directed only on the 
passage of water resistance. 
Controls in the Joints are com 
parable in value with 
the control in a terrestrial 
* apparatus. Finally, the amplitude 
of oscillations of the body i s 
not great (~ 3°). 
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The noted difference In the dynamics of an underwater and terres- 
trial apparatus does not always take place. The conditions In which 
the dynamic characteristics of the underwater apparatus bear a signif- 
icantly oscillatory nature, and also the mechanisms of these oscillations 
will be described below. 

As calculations show, the decrease of the distance between 
the point of support and the center of masses leads to an increase of 
the frequency of oscillations of the dynamic characteristics: In that 
temporary section appear a great quantity of mlnimums and maxlmums 
of functions. This tendency of Increase of the frequency of oscill- 
ations with a decrease of p is seen In the picture of Pig. 8. In 
this figure, results are shown of calculation for a very small value 

-C 

p ® 10 ^m (al] remaining parameters are as In the base variant). It 
Is possible to explain the noted effect theoretically also (In a 
linear approximation) . From the formulas for the frequency in the 
natural oscillatlono of the body wy in (5.2) we find: 

where C^, are positive. Independent of p^, constants. From here 
it follows that with a decrease of p the frequency of natural oscilla- 
tions of the body of the apparatus increases, and this means that the 
frequency of oscillations of the dynamic characteristics caused by it 
also Increase. 

In Fig. 9 the results are shown of calculation of a variant, 
the Initial parameters of which are distinguished from the parameters 
of the base variant only by the velocity of motion: Vjj = 0.2 m/sec. 

The change of velocity of motion acts little on the frequency of 
natural oscillations, which is evident from comparison of Fig. 9 with 
Fig. 7: for one and the same interval an identical number of extremes 
of functions take place. Increase of the number of oscillations in 
a period of the length of time of a step is explained by the increase 
of this very period. /4o 
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Figure 8 . The dynamic characteristics 
of the apparatus in its motion on one 
step. = 3180 kg, Pg = 3130 kg, 

Vn “ 0.5 m/sec, R .= 0.9 m, = 0.6 m, 
** u C 9 o 

= 0.0^1 m. 


The amplitude of 
oscillations of the body some- 
what increases, but remains 
small: the initial value Yq “ 
- 0 °. 6 , and the amplitude 
of oscillations near this 
value have an order '•^O. 

With an Increase of the 
weight of the apparatus, the 
maximum values of the forces 
of reaction and the amplitude 
of oscillations of the body 
arise, which Pig. 10 illus- 
trates (this supports also 
the qualitative analysis of 
the motion of the apparatus 
carried out in the boundaries 
of linear theory). In this 
figure the initial parameters 



r — 

Figure 9 . Dynamic characteristics of the apparatus in its motion on one leg. 
P^ = 3180 kg, Pg = 3130 kg, = 0.2 m/sec, » 0.9 m, p^ » 0 . 6 m, p^ = 0 . 3 m. 
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are dlstingulahed from the parameters of the base variant only by 
the weight of the apparatus ■ 3330 kg, suoh that there Is a 
difference of weight and buoying force P^-Pg ■ 200 kg. The maximum 
value of the horizontal component of reaction here becomes ~ 100 kg, 
the maximum amplitudes of oscillations of the body ■<' 20°. 

The oscillatory nature of the dynamic characteristics Is Induced 
by the fact that In the examined case 

and the summated momentum of the force of gravity and Archimedean 
force tends to put the body of the apparatus In the position y* 0 . 
Oscillation of all the charcterlstlcs Is a reflection of the natural 
oscillations of the body near the position y “ 0. 



Figure 10. Dynamic characteristics 
of the apparatus In Its motion In 
one step. ■ 3330 kg, Pg ® 3130 

kg, Vjj ■ 0.5 m/sec, » 0.9 m, 

p ■ 0.6 m, p ■ 0.3 m. 

o c 


It Is possible to try 
to attain a planar, almost 
nonosclllatory change of 
reactions and controls In time 
T of the length of one step, 
having selected this length 
sufficiently small In comparison 
with the period Tc of the 
natural oscillations of the 
system. So, for example, 
having Imposed the condition 

T<±Tc , (6.2) 

It is possible to expect that 
the dynamic characteristics In 
time T will not be able to 
complete a great number of 
oscillations and will not have 
more than one or two extremes. 
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In terms of the parameters 
of control of small oscillations 
(5.1) the condition (6.2) 
is recorded in the following 
form: 

The characteristics of 
motion of the apparatus with 
a velocity Vjj ■ 0.97 m/sec, 
which corresponds to the sign 
of the equality in (6.3) (with 
values of the remaining para- 
meters of the base variant) is 
shown in Fig. 11. It is obvious 
that the dynamic characteristics 
are, in this case, actually 
weakly oscillating. 


Figure 11. Dynamic characteristics 
of the apparatus in its motion in 
step. = 3180 kg, Pg = 3130 Kg, 

Vji = 0.97 m/sec, 0.9 m. 


= 0.6 m. 


0.3 m. 


The variants examined 
above answer the case of (6.1). 
If the parameters are selected 
in such a manner that (6.1) is 
not fulfilled, then the picture 
of change of the function 


characteristics will be such as in Fig. 12. 


0.6 m, » 0 . 5995 m) . 


"^o ’ c 

Here, as in the base variant, at the moment of change of the step, 
the point of support projects on the support point. It is obvious 
that the characteristics change very smoothly in contrast to the case 
of fulfilling the conditions of (6.1). Let us note, however, in the 
examined case, the solution found is unstable, as in the tasks of 
terestrlal walking. (Instability may have a positive consequence for 
the controllability of the apparatus.) 


As has been shown earlier, at least two periodic modes of 
operation take place: in the area y “ 0 and in the area y * it . The 
latter mode for those values of parameters, which are also in Fig. 12 










Figure 12. Dynamic characteristics 
of the apparatus in its motion on 


one step. 


3180 kg, P„- 3130 


kg, Vjj- 0.5 m/sec, ■ 0.9 m, 
p_«0.6m, p_" 0.5995 m. 


Figure 13 . Dynamic characteristics 
of the apparatus in its motion 
on one step. 3180 kg, 

Pg 3130 kg, Vu ■ 0.5 m/sec. 

R„. ■ 0.9 m, p^ ■ 0.6 m, 

P ■ 0.5995 m (solution in the 

c , 

area y ■ , y ■ 0) 


area y ■ , y ■ 0; 

are shown in Fig. 13 . Let us suggest that stability in the area of 
Y ■ 180® depends on the sign of the value 

stability in the area y “ 0 only on the sign of the value P_p -P.P . 
Therefore, instability in the area of y ■ 0 does not certainly answer 
the stability in the area y ■ 180°. Both modes may appear unstable 
(such as takes place in the cases shown in Fig. 12, 13). 


In Fig. 1^1 the case of a stable mode is shown in the area y ■ 180°. 
The parameters are taken from the base variant, which give a stable 
motion in the area y ■ 0. Stability is kept also in the area y ■ 180° 
owing to the fact that J1 ■ J„-Mpn^H<0, as this was explained in the 




Figure 1^1. Dynamic characteristics 
of the apparatus in its motion on 

one step. ■ 3180 kg, 3130 kg, 

Vjj* 0.5 m/sec, 0.9m, 0.6 m, 

■ 0.3 m (solution In the area 

Y ■ IT, -Y « 0) 


preceding paragraph. Very strong 
oscillations of all character- 
istics are noted — the 
consequences of a small (by the 
module effective, "Introduced”) 
momentum of Inertia Jjj ) . 

It Is useful to recall that 
stability (Instability) here Is 
understood In the sense In- 
dicated above, in Section 5» 

— as stability of the angular 
motion of the body with in- 
variability of the forward 
motion and motion of the legs. 

In conclusion It Is neces- 
sary to note that the disregard 
of the combined masses qualita- 
tively changes nothing (all rules 
remain in force). Only quantita- 
tive relations are able to 
change. What has been said 
directly follows from the 
equation of motion of the 

where the 


apparatus. This Is obvious from comparison of Fig. 7 and 15 
results of calculations are shown of one and the same variant with 
regard (Fig. 7) and without regard (Fig. 15) of combined masses. 

Conclusion 


The main quantitative results of the calculations done are 

given In two consecutive tables. The maximum values, rounded off, 
of the values of reaction and controls are combined in them, and also 
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Figure 15* Dynamic characteristics 
of the apparatus In its motion on 

one foot. • 3180 kg, Pg ■ 3130 

kg, Vg ■ 0.5 m/sec, • 0.9 m, 

■ 0.6 m, p ■ 0.3 m. (Solution 
o c 

without calculation of combined 
masses) 

TABLE 1 


the median values and amplitude 
of oscillation of the angle of 
Inclination of the body of the 
apparatus. Shown alec is the 
stable or unstable angular 
motion of the body In that 
limited sense which was ex- 
amined In the work (with an 
Invariable forward motion). 

The last graph of the table 
contains the conditional number 
"n" characterizing the smoothness 
of the dynamics of the apparatus: 
this Is the "number of peaks" 
(maxlmums and mlnlmums) In the 
controls and the horizontal 
component In the time of the 
length of one step. 

Table 1 contains the In- 
dicated values for the apparatus 
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with an effective weight ■ 50 kg. with various values of velocity /3.1 
of the apparatus. In Table 2 with a fixed value of velocity (Vjt ■0.5 
m/sec) the effective weight of the apparatus Is varied. 


TABLE 2 

Taoiica 2 








mm 

50 

SX yes 

50 

30 

40 

0®i 3“* 

. 1 

200 

a, no 

220 

ZOO 

250 

00*20® 

_j_J 


It seems natural to consider such motion "good", which, with 
the other equal conditions provides a) a minimum of maximum values 
f't the dynamic characteristics and b) a smoothness of change with 
time of these dynamic characteristics. Actually fulfullment of con- 
dition a) Is advantageous energetically, and fulfillment of condition 
b) represents a simplicity of control by the apparatus In a nominal 
mode of operation. In aggregate fulfillment of conditions a) and b) 
make It possible to call it "good controllability" of the apparatus. 

From the given tables it is obvious that the motions with a 
moderate speed meet these conditions, with an upper position of the 
center of masses relative to the point of support of the legs and 
moreover unstable in the sense introduced above. Stability mars the 
power engineering, particularly with a lower (y*180®) position of 
the center of masses. In addition stability leads to a strong 
oscillation of the dynamic characteristics, whereas the unstable modes 
give smooth, monotonic functions. 

Instability of the angular motion In the sense examined here 
does not present any danger. It is known that adequate changes of 
the law of control are able to stabilize similar instabilities. Stabili- 
zation is achieved for example, by a small systematic change of the 



















length and length of time of the step [7] 


Prom the second table it Is also obvious that an Increase of 
apparatus weight, naturally, leads to a corresponding Increase of 
maximum values of the dynamic characteristics. 

The conducted research of a model task of the dynamics of /53 

underwater two-legged walking allows the following conclusions to be 
made. 


1. Modes of operation of motion exist which provide Its 
conformablllty with periodic compensating oscillations of the body 
of the apparatus. 

2. At least two classes of periodic oscillations of the body 
exist: with an upper and lower position of the center of masses of 
the body relative to the point of support of the legs. 

3. The conditional stability of the angular motion of the body 
Is understood In Invariable forward motion of the point of support 

of the legs. If the maximum momentum of the Archimedean buoying force 
Is greater than the maximum force of gravity, then the reducing momen- 
tum Is positive and the "upper" mode of motion In the Indicated sense 
Is asymptotically stable. In an opposite case — unstable. 

Stability or Instability of the "lower" mode ir defined not 
only by the sign of the reducing momentum, but also by the sign of the 
Introduced momentum of Inertia. The asymptotic stability Is achieved 
when these values have different signs. 

5. The oscillatory nature of change of the dynamic characteristics 
of walking is accompanied by the stable mode of motion. These 
oscillations take place with the frequency of the natural oscillations 
of the apparatus. A smooth (often monotonlc) nature of change of the 
dynamic characteristics is accompanied by an unstable mode of motion. 







6. With the change of legs the values of the support reaction 
and of controls undergo a disruption which Is a consequence of single- 
support movement . 

7 . The vertical component of the support reaction at the point 

of support Is little distinguished from the difference between the 
force of the weight and the buoying force of the body. The horizontal 
component of the support reaction has this order: It Is less than 

the vertical component only by one and one half to two times. Almost 
henceforth the horl.o'.tal component of reaction Is positive, there /52 
being a strong qualltat !■'. ,? and quantitative difference in the behavior 
of the horizontal component of the underwater and equivalent terres- 
trial apparatus. This is explained by the fact that in underwater 
walking the horizontal component of reaction must compensate for the 
resistance of the envlornment. 

8. The controlling momentums In the Joints of the legs in order 
of value are such as the controlling momentums for an equivalent 
terrestrial apparatus. In the support leg the momentums in the hip 
Joint are approximately one factor of that of the momentums in the 
shank. The momentums in the moving leg are small in comparison with 
the momentums in the support leg. In addition, the momentum in the 
hip Joint on one leg, as a rule, changes sign while the momentum 

in the knee Joint in_ basis of sign does not change sign and is positive. 

9 . The unstable modes require a lesser maximum value of control- 
ling momentums in comparison with stable. 

10. The unstable mode with an upper position of the center of 
masses of the body by the minimum of maximum values of the dynamic 
characteristics and by the smoothness of their behavior is most 
preferable. The worst in this sense is the stable mode with a lower 
position of the center of masses of the body. 








Appendix; Generalized Method of Cords 


m 


In the present work, for the solution of the boundary task 
(1.^4)-(1.6) , the numerical, Iterative, generalized method of cords 
was used. The basic Idea of this method Is Included In the follow- 
ing. Let It be necessary for us to find the solution of a system 
of transcendental equations 

/fx) • 0 , 

(PI) 


where f (x) * f^ (x) , . . . ,f^(x) Is the vector-function of the vector 


argument x ^ =(x^ , . . . ,x^) of the dimension n (here and henceforth 


by the prime Is meant the operation of transposition). Let us assume 


,(n) 


of the solution of the system 


that we know several approximation x 

^ let us approximate the nonlinear state- 


(Pl). In the population x 
ment f(x) by the linear 

/(jc) a A(j:- X *■ a 


(P2) 


Here A Is (nxn) — the matrix, and B Is n — the vector. Let us find 
the unknown matrix A and the unknown vector B from the following 
conditions . 


Let us assign In the population x^’^^ the point x^^^(l = 0 , 1 , . . . ,n)and 


compute the value f(x^^^) (1 = 0,1, . ,n-l) . Then let us require that 

( 1 ) 


at the points x (1 = 0,1,..., n) the values of the nonlinear and 
linear statements coincided 


f(x^'*) - A(x , i»0, /, .... n 


.tit tnt 


This Is equivalent to the matrix equations 

Ant 


where AF and Ax are (nxn) — the matrixes of the form 


(P3) 


aP = 



Af 


‘fr 

■‘j£ • • ' 
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A JE^ 4JC|g 


axM 


/5^ 


£>Xn • • • ^«f/t 


^V •v^'* . 

Let us assume that the points are selected such that 

det ^ 0. 

Then from (P3) we find 
A ■ AP(AX)‘ 

which together with (P2) and (P3) gives 


.-1 


f(x) * AF(AX)"^(x-x^"^) + f(x^”h 


The following approximation is received in the usual manner 

by means of solution of the system of linear equations 


f(x •”) - o 

Supposing that 

det AP ^ 0 

we find 


Having received let us compute fCx^*^^^^). Then among 

the points x^®] x^^\ . . . let us eliminate the most remote from 

jj(n+l) pQj^ptg^ ipj^g distance between points is considered in several 
assigned metrics. In the given work the distance p(x^^^,x^^^) between 
points x^^^ and x^^^ were computed by the formula 


fifx'.a:' 






where Pj are the assigned nonnegative weight coefficients, 
of the eliminated point 1^ is defined by the condition 


y,rx“'ix"^} 


• max plx”\ x""'j 

Oikitn X * ' 


Tho number 


The following approximation of x^'^^^^is calculated by the same formulas 
that were used for but only in these formulas it is necessary 

to substitute x^^*^for x^*^^ (concurrently f(x^^*^) for f(x^*^^), if 
Iji ^ n, x^*^^ for (respectively f(x^*^ ) for f(y'^'*^^ ). Having /55 

obtained in an analogous manner let us calculate then 

x(n+4) forth. The iterative process ends if at the s iteration 
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the conditions are fulfulled 


* l*Sj ; s^ft » 

where are the assigned positive numbers which characterize the 
precision with which it is necessary to know the solution. 


The weight coefficients Pj are selected from the considerations 
of dimensions. In particular it is possible to set 

4 ••• f 

if Xj is dissimilar, or 
if Xj is similar. 


It is possible to show that the generalized method of cords 
converges if the zero approximation (points x^®\ x^^\..., x^^^ 
are sufficiently close to the solution. Here the speed of convergence 
of the method in proportion to the approximation to the solution 
strives to the speed of convergence of the method of Newton. Practical 
use of the generalized method of cords showed that the general volume 
of computation which must be done in order to receive a solution 
with the assigned precision is less than with the use of the method of 
Newton, despite the fact that the number of Iterations which must be 
done is larger. The generalized method of cords by speed of convergence 
being no worse than the method of Newton, has one significant advantage: 
with its use it is not necessary to compute the matrix of partial 
products - » 

9x 


which for many tasks has a considerable value. In addition, in contrast 

to the method of Newton, when for each iteration it is necessary to 

convert the matrix of n factor djL , the iteration process in the 

dx 

generalized method of cords can be so limited, that on the first 
iteration the matrix will be converted by an n factor, and in all the 
following — the matrixes only of the second order, which is very /5( 

H3 




I ini w i i wn 


” "rf • T»TT^ » * V 




Important If n is Important. The latter Is ach^e,ved in the following 
manner. Let us designate through the matrix AP. which Is 

computed In the first Iteration by the points x^®\ x'^\..., 
see formula (P4). In the first Iteration In the process of obtaining 
the first approximation It is necessary for us to convert the 

matrix of n factor AP^^^. According to the generalized method of 
cords, after the first approximation of x^”'*’^^ls obtained# one of 
the points x^^\ x^^\..,, Is eliminated. Let the eliminated 

point prove to be the point with the number 1. 'Then In the following 
iteration it will be necessary for us to convert the matrix 


^ A • • • i • • • ^ 


where 




is the J column of the matrix AP 
In the form 


( 2 ) 


The matrix AP 


( 2 ) 


may be recorded 


Here 

s'*‘\ 



^cn*o ^ 



foo . 

. . 010 

...OK 



(// ■ . 

. /// 

... tJ 


For the conversion of matrix AP 


( 2 ) 


we use the formula 


(A*&J)C)‘'^A'-A'b(D* CASf'cA*. 


(P5) 


This formula takes place for any matrix A, B, C, D for which all 
operations are fulfilled, which are figured In (P5) (it is easy to 
be assured of the accuracy of (P5) by direct checking). Supposing 
In (P5) 

c.c"] 


we obtain 




AP 


( 2 ) 


Prom this formula It is obvious that conversion of the matrix 
is reduced to conversion of the matrix of the second order 
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E2 + In as much as the matrix AP^^^”^has already 

been computed In the preceding iteration. In an analogous manner in 
each final Iteration, using the results of the preceding iteration, 
the conversion of matrix AP will reduce to conversion of the matrix 
of the second order. 

In the particular case when i= n, that is when point is 

( 2 ) 

eliminated, conversion of the matrix AP may be reduced to division 
by the number. Actually in this case 

n 

from this it follows 

T^CAPf**'' B 

Thus when the point which answers the approximation obtained in the 
pi-eceding iteration does not take part in the construction of the 
given approximation conversion of the matrix AP reduces to multipli- 
cation of the knowns of the matrix and division by the number. 
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